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Abstract 

A Fermi's Golden Rule for population transfer between instantaneous eigenstates of elliptical 
quantum billiards with oscillating boundaries is derived. Thereby, both the occurrence of the 
recently observed resonant population transfer between instantaneous eigenstates [F. Lenz et al. 
New J. Phys., 13, 103019, 2011] and the empirical criterion stating that these transitions occur 
when the driving frequency matches the mean difference of the latter are explained. As a second 
main result a criterion judging which resonances are resolvable in a corresponding experiment of 
certain duration is provided. Our analysis is complemented by numerical simulations for three 
different driving laws. The corresponding resonance spectra are in agreement with the predictions 
of both criteria. 
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I. INTRODUCTION 



Classical driven billiards of varying geometry have been subject to intensive research dur- 
ing the last years []]-13|. While billiards are, in general, important models to study aspects 
of nonlinear dynamics, semi-classics or (quantum) chaos 14|, ll5j, driven billiards addition- 
ally facilitate the study of non-equilibrium dynamics. As one of the key topics concerning 
driven billiards, Fermi acceleration (FA) and the related conditions for its occurrence have 



gained much attention js-ll3|. FA describes the unbounded growth of energy of particles 
that repeatedly interact with a time-dependent potential that is usually modeled by a mov- 
ing billiard boundary and was originally proposed by Enrico Fermi as a possible mechanism 
to explain high-energetic cosmic radiation [16|. The infamous Fermi- Ulam model (FUM) is 
basically a one-dimensional billiard with a moving boundary and it was found that FA is 
present in the FUM only for non-smooth driving laws 13] . The general conditions for the 
occurrence of FA are still under debate. Originally, it was assumed that a sufficient condi- 
tion for the occurrence of FA in a driven two-dimensional billiard is the presence of chaotic 



131 ] . However, it turned out, 



regions in the phase space of the corresponding static billiard 
that driving an ovally shaped billiard which has a mixed phase space in a certain mode 
does not lead to FA 18]. On the other hand, it was shown that the classical driven ellip- 
tical billiard does show FA although its static counterpart is completely integrable 0, 10]. 
Furthermore, while correlated motion suppresses FA for smooth driving in the FUM, it was 
found that correlations can even cause exponential FA for smooth driving laws in a related 
two-dimensional model 
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121. 



Although it is known that periodically driven quantum billiards with a discrete Floquet 
spectrum can not exhibit FA [191 ]. it is natural to complement the study of the classical 
dynamics of a system by analyzing its quantum behavior. While one finds many studies 
to the quantum version of the one- dimensional FUM (see 2CH22I] and references therein), 



literature is very sparse on driven quantum billiards of higher dimensions 231-125]. 



In particular, [23J presents a method to solve the time-dependent elliptical quantum 
billiard. The main result was the numerical observation of resonances in the population 
transfer probability between instantaneous energy eigenstates. These transitions could be 
reproduced in an effective Rabi-model and captured by a criterion stating that resonances 
occur whenever the difference of corresponding time-averaged energy eigenvalues matches 



an integer multiple of the driving frequency. However, an explanation for this criterion was 



not given in 



23]. 



Here, we develop a systematic perturbative analysis of population transfer for the system 



analyzed in 



231 ] and a generalized driving law. In this framework a Fermi's Golden Rule 



26| is derived for elliptical quantum billiards with oscillating boundaries which explains 



the key observations in [23], i.e. the occurrence of resonant population transfer between 
instantaneous eigenstates and the empirical criterion relating these resonances with the 
spectrum of instantaneous eigenstates and the driving frequency. As a second major result, 
we provide a criterion to decide whether a predicted resonance can be resolved in a possible 
'experiment' of a certain duration. Finally, the numerical studies in [23} are complemented 
by a corresponding analysis of further driving laws. The predictions derived within our 
perturbative analysis will be shown to provide a perfect agreement with the numerical results 
in all cases. 

This work is structured as follows: Chapter III Al provides a short summary of the so- 
lution of the time-dependent elliptical quantum billiard as developed in 23j, followed by 
transformations that bring the Schrodinger equation into a form being convenient for the 
application of time-dependent perturbation theory. In chapter III we will calculate the 
transition rate between two instantaneous eigenstates per unit time is calculated in first 
order perturbation theory and an approximate population dynamics in the near-resonant 
case is derived. We find a criterion for the resolvability of predicted resonances in a possible 
experiment of certain duration. Finally, in chapter II I II we present and analyze numerical 
results for three different periodic driving laws. 



II. TIME-DEPENDENT ELLIPTICAL BILLIARD AND ITS ANALYTIC TREAT- 
MENT 

In the following, we first summarize our approach to a numerical solution of the time- 



dependent Schrdinger equation of the elliptical billiard as presented in [23}. We will trans- 
form the time-dependent Schrodinger equation (TDSE) into a convenient form and finally 
develop a perturbative approach for the periodically driven billiard in the second part of 
this chapter. 
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A. Setup 



A wave function ty(x,t) in a driven elliptical billiard obeying the TDSE 

h 2 

ihd t V(x,t) = -— Att(£,f), (1) 

is subject to Dirichlet boundary conditions, ^(x,t)\gB = 0, on a boundary OB of elliptical 
shape: 

dB = 1.x = (x,y) J e ) 

Here, the semi-axes of the elliptical boundary, a and b, are assumed to be arbitrary smooth 
functions of time, i.e. a = a(t) and b = b(t). 

The time-dependent boundary conditions can be handled by a coordinate transformation 



23] 



-r) • ) W i <3) 



y I V V V W) I \ y 

that maps the time-dependent elliptical boundary onto a static boundary of the shape of a 
unit circle. Applying (EJ) to together with a unitary transformation, 



and extracting a volume-dependent prefactor a(t)b(t) from the wave function ^(x, t), we 
are led to an effective SE 

ihd t A( V ,Z,t) = H e ( V ,Z,t)A(r ] ,Z,t), (5) 



where the effective Hamiltonian H e contains time derivatives of the prefactor a(i)b(i) and 
of the unitary transformation U of the left-hand side of the TDSE: 

£ - t] ^ if iw) w- + wM + 5" ( a(t)a(i) " 2 + b{tme ) ■ (6) 

The introduction of the unitary transformation (jlj) ensures that H e is Hermitian. Due to 



the extracted prefactor a/ a(t)b(t), the effective wave function 

A{V,£,t) := '(^O,*)*^" '(^O.*) (7) 

is normalized to 1 on the domain boundary of the unit circle C := {x = (x, y) J G M 2 \x 2 + y 2 < 1 } 
and the coordinate transformation ([3]) makes A subject to the Dirichlet boundary condition 



ty(x, t)\gc — 0. The reader is referred to 23[] for a similar, more detailed derivation of the 



effective Hamiltonian and equations of motion. 

A complete set of orthonormal functions on C is given by the eigenfunctions of the static 



circular billiard 
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3l|: 



$ n , m ( P , 0) = 1 r 4(fcm,nP) eM - ( 8 ) 

V J m+l {km,n ) 

p and (f) can be calculated from rj and £ by rj — p cos and £ = p sin 0. J m is the cylindrical 
Bessel function of order m and k m>n is its n-th root, n and m are called radial, resp. angular, 
quantum number for obvious reasons. If we expand the effective wave function A in terms 
of the eigenfunctions of the static circular billiard, the effective SE (J5]) becomes a linear 
homogeneous ordinary differential equation of first order in time and can, thus, be solved 



numerically by standard methods 23 1. 



A main result of 23] was the observation of resonant population transfer between so-called 



instantaneous eigenstates of 

-h 2 f l d 2 l d 2 \ 

M ~ 2p \a 2 (t)d V 2 + b 2 (t)de) ' 

We understand instantaneous eigenstates as follows: The semi-axes a and b are parameters 
of Hm that change in time. If we evolve our system solely by Hm in the SE, start the 
system in an initial state that corresponds to an eigenstate of Hm at t = and change a 
and b sufficiently slowly, then we define the instantaneous eigenstate of H M at time t as the 
time-evolved wave function of the system at time t in accordance to the adiabatic theorem 



of quantum mechanics 



29] 



The Hamiltonian Hm (EJ) is part of the effective Hamiltonian H e ([6]). Its complementary 
part is 

Hp = H e — H M = \p (a(t)ti(t) v 2 + b(t)b(t)£ 2 } . (10) 

Population transfer between instantaneous eigenstates of Hm takes place by two different 
mechanisms in the billiard. First, as a and b are of course not changed sufficiently slowly, 
diabatic population transfer between the instantaneous eigenstates of Hm will take place. 
Additionally, the Hamiltonian Hp triggers population transfer as it is non-diagonal in the 
basis set of instantaneous eigenstates of Hm- 

Introducing the volume of the elliptical billiard, V(t) = a(t)b(t), and the ratio of the 



semi-axes, r(t) = b(t)/a(t), Hm can be rewritten in the much more convenient form 



1 / d 2 Id 2 
M{r):=--(r— + -—) (12) 



where we call 

1 f d 

M{r) := -- (^r- 2 , ^ Q ^ 

the Mathieu operator as its eigenfunctions are just ordinary and modified Mathieu functions 
as they appear in the solutions of the static elliptical billiard. If we label the eigenstates of 
M(r) by |n; r) with eigenvalue g n ( r )> \n; r(t)) are, of course, the instantaneous eigenstates of 
Hm and E n {t) = j^m Qn{ r {t)) t ne corresponding instantaneous eigenvalues of H M , i-e. we 
have 

H M = ^n;r(t))^M(n;r(t)\. (13) 

Note that M(r) is invariant upon sign-change of r] and £. One can therefore choose its 
eigenstates \n\ r) such that they are also eigenstates of the parity operators that change the 
sign of rj or £. In this context, we will refer to \n; r) having even or odd rj-, resp. parity. 
Note that the effective Hamiltonian H e is also invariant upon sign-change of rj and £ and, 
consequently, only couples instantaneous eigenstates that have the same rj- and £- parity. 
The Hilbert space therefore splits into four uncoupled Hilbert subspaces. 
We choose the following ansatz for the effective wave function A, 

|A(t)) = ^c n (t)e-^W|n;r(t)), (14) 

n 

with time-dependent expansion coefficients c n (t) and 

Mt) ■= j dt ' E n(t') = f dt '^) ft»( r (0) ( 15 ) 

being the time- integrated instantaneous eigenvalues of H M - If we put this ansatz into the 
the SE OS]) and note that d t |n; r(t)) = f(d r \n; r))\ r=r (t) = rd r \n; r(t)) , we get a SE for the 
coefficients c n (t): 

ihd t c n (t) = c m (t)e"s^W-^W) ((n; r(t) \ H F (t) \ m; r(t)> - ihr(t) (n; r(t) | 9 r | m; r(t)» 

m 

(16) 
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For non- degenerate eigenstates \n; r) and \m; r), one can express the second matrix element 
on the right-hand side of ffTB"]) as 

' (n;r\(drM(r)nrn-,r) for ^ 

(n;r|«9 r |m;r> = <[ r (17) 

(n; r | (9 r | n; r) for n = m. 

It is now interesting to notice that the representation of M(r) in the eigenbasis of the static 
circular billiard (JBJ), (3> n ' )m > | M(r) | $ n ,m), is not only a Hermitian, but a real symmetric 
matrix (cf. Eq. ()A14j) in appendix E| . We can therefore choose the expansion coefficients 
of the eigenstates of M(r) in the eigenbasis of the static circular billiard, {<& n ',m' \n;r), to 
be real. It follows that also the expansion coefficients of d r |n; r), (3> n / OT / | 9 r | n; r), are real. 
Thus, 

(n; r | <9 r | n; r) = ^ (n; r | $ n /, TO /) ($ n ',m' I I ™5 r) (18) 

n',m' 

is also real. On the other hand, due to normalization of the eigenstates |n;r), expression 
(|18p has to be purely imaginary and is therefore identical zero. Noting that 

d 2 1 d 2 



d r M(r) 



drj 2 r 2 d£ 2 
1 / d 2 Id 2 



r- 



r \ dr] 2 r d£ 2 
M(ir) 



(19) 



ir 

further simplifies Eq. ( JT7|) . From now on, we will restrict ourselves to periodic driving laws, 
i.e. a(t + — ) = a(t) and b(t + — ) = b(t), and cases where all populated instantaneous 
eigenstates are non-degenerate. These conditions also include all cases that are discussed in 



231 ] . We further introduce a rescaled dimensionless time r := j-t and finally put Eqs. ((T 



and (USD back into ([161). 



^ T c n (r) = J]c»We"^ ( *" (rH " (T,) ^ (n;r(r) |ffr(T) |m;r(r)> 

E, , a-™/,,, j, r(r) (n; r(r) I M(zr(r)) I m; r(r)) , , 

c m r e--(^^-^w) v ' 1 ;' 1 ' 1 ;/ . 20 
m ^ n r ( r ) q n (nr)) -q m (r(r)) 

We point out that the modulus of the first term in (1201) depends linearly on the driving 
frequency u> while the modulus of the second term is independent of u. We therefore ex- 
pect the first term to be dominating for large driving frequencies, while the second one 



should be dominant for small driving frequencies and should especially couple neighboring 
instantaneous eigenstates due to the denominator g n (r(r)) — q m {r(j)). 

Obviously, due to periodic driving, all terms on the right-hand side of ( l2"0l but the coeffi- 
cients c m (r) and the phase factors exp(— 27ri(0 m (r) — n (r)) / Jtw) are one-periodic functions 
in r. It is therefore possible to represent them by discrete Fourier transforms. Before we 
do so, we split m (r) — n (r) into a non-periodic part hu mn r := (0 m (l) — n (l))r and a 
one-periodic part /j,Az/ m „(r): 

0m (7) - (t>n{r) = hv mn T + hAu mn (r). (21) 

We then combine the one-periodic phase factor exp(— 2^iAv mn {r) / uo) with the other one- 
periodic terms on the right-hand side of (120]) and Fourier transform the results: 



1=00 

u J2 if m e- 2rfr = e ~^ AiWT) ^ ( n ; r ( r ) I h f{t) I m; r(r)) (22) 



i=— 00 



27TI , 



JM^ (^)iy))Mr)) forn^m 
2=-oo I for n = m. 

Before we put ( 122|) and ( )23|) back into the SE ( 120|) . it is useful to perform a unitary 
transformation, c n (r) = exp(— iuFQ n T)b n (r). Note that we do not include D™ in the unitary 
transformation as it is zero by definition ( 1231) and that Fq u is completely independent of u. 
This unitary transformation, together with Eqs. ( 1221) and ( 1231) . leads via ( 1201) to a SE for 
the coefficients 6 n : 

*6«(t)= E e^ ie r-{uFr + Dr)bm{r), (24) 

m^n for i=0 

where we have defined the abbreviation 

f:= — + — (F nn -F mm )-/. (25) 

OJ 27T 



Note that the solution of Eq. ( 1241) determines the complete physics of periodically driven 
elliptical quantum billiards. 

B. Perturbative analysis 

We will now use time- dependent perturbation theory (TDPT) to find an approximate 
solution of Eq. fl24|) in first order. To do so, we formally affix a parameter A to F™" 1 and 



S 



D™ m to keep track of the order of perturbation and will set A to 1 at the end of our calculation: 

F nm = X . pmm Jjnm = X . pnm An expansion f ^( r ) in \ gives: fr n ( r ) = J2™=0 ^^(r) . 

As A should track the order of perturbation, it is natural to choose the initial values of bn^ 
as bn\o) = <5 Pi o " b n (0). Inserting this ansatz into (12^1) and equating equal powers of A yields 
up to first order: 

ifoW = = const. = b n (0) (26) 



ifeW= J] e M " mT K m + Dr)6 m (0) (27) 



mj^ii for Z=0 



6 ? ) W= E O -flnm (^ m + iJP m )M0) ■ (28) 



m,l 
m^n for 1=0 



1. Population transfer probability 



We are now able to calculate the population transfer probability between two instanta- 
neous eigenstates which will lead to a systematical understanding of resonant population 
transfer as it was observed in {23)]. For this purpose we assume that the wave function 
I A) was initially in the (undriven) eigenstate \k;r) and then calculate the evolution of the 

population of the eigenstate \n; r) (n 7^ k). Population transfer in first order gives: 

27Tie? k r _ 1 p-2-Kiefr _ 1 
:= |6«(r)| 2 = £ ——{u,F«* + D*){ U F?* + D* v nk ). (29) 

We would like to calculate a population transition rate per unit time from ( 1291) which is de- 
fined as r^ fc := lim^oo p l n fc (r)/r. Note that (e 27ri0 "^ - l)/2n&? k = ie^ T sin(7r0p*r)/(7r0p*) 
grows linearly with r for 8f k = while it oscillates periodically with an amplitude 1/7T0" 
(which is independent of r) for 7^ 0. Due to 9f /k — 9™ k — I' — I, we can therefore neglect 
all terms in (|29|) with I 7^ /' for r being sufficiently large. 



1 1 



nk 



lim 



Pnfc( r ) 



D 



nk 1 rv'<'l" 

nfc I ^ 



(30) 
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By applying appropriate transformations, we have handled the time-dependent boundary 
conditions of the billiard by introduction of a time- dependent external potential. This 
enabled us to derive Eq. ( 13 Op which is a Fermi's Golden Rule 26] for driven elliptical quantum 
billiards. It states that efficient population transfer in first order between the instantaneous 
eigenstates \k;r) and \n;r) is only possible for = 0. We can now use (1251) to calculate 
corresponding resonance frequencies, 



hl _ I ± y/Z 2 - Av nk 5F£ k 

"res ~ 2«F? ' (31) 

where 2tt ■ 5FQ k := Fq u — F fcA: has been defined. Numerical experience shows that 5FQ k 
is usually a very small quantity. The "+"-term in ([3~Tj) thus corresponds to a very large 
resonance frequency. Restricting ourselves to not too strongly driven billiards, we will neglect 
this term from now on. If we develop the "—"-term in (l3Tj) about 5FQ k m and use the 
definition of v nk in Eq. ( 12 ip above, we find: 



I ■ ujfj = v nk = f E n {r') - E k {r')dT' 
Jo 



(32) 



Thus, only when the one-period average difference of two instantaneous energy eigenvalues 
matches an integer multiple of the driving frequency, resonant population transfer between 
the corresponding instantaneous eigenstates can occur. This is precisely the empirically 



found criterion in 



23| and has herewith a theoretical basis. The result justifies to call the 



Fourier summation index / "photon process order" of a population transfer in analogy to 
the interaction of light and matter. 



2. Applicability of first order TDPT 

Not all predicted resonance frequencies ( |3T|) are of equal importance with respect to their 
experimental observation and we will now derive a criterion to discriminate them. In the 
resonant case 9f" k = 9^ k = 0, Eq. ( |29|) reduces to: 

Pn k = T 2 Y,\" F i nk + D ? k \ 2 ( 33 ) 
I 

The reader is reminded that (J33J) only holds for n ^ k, while for n = k, p\ k = holds 
as 6 1 = implies / = and this case just had been excluded from the summation in 
Eq. (T28|) . Consequently, the instantaneous eigenstate \k;r) gets exclusively depopulated in 
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first order TDPT. We can therefore calculate the time T- mt at which the population pk of the 
instantaneous eigenstate \k;r) gets negative and therefore unphysical: 

PfcOint) = 1 - Pnk( T int) = (34) 



'int 



yak 



(35) 



T^oK"* + D i 

The summation index 9f k = in (1341) and (135*]) means that it should only be summed over 
states n and photon process orders / that satisfy the resonance condition 9f h = 0. This 
means in all practical examples that the sum only consists of a single term. 

ri n t is a measure for how fast a population probability transfer takes place. It is, thus, 
reasonable that we will not be able to fully resolve resonances that correspond to an inter- 
action time T int that is much larger than the actual runtime r mn of a possible experiment. 
In this case, population transfer will have been stopped before the maximal theoretically 
possible amount of population probability will have been transferred from one instantaneous 
eigenstate to the other and our ability to resolve a resonance in corresponding observations 
is diminished. 

On the other hand, we understand that the transition rate (1301) has been calculated in 
the limit r — > oo and the included 5-function is the result of a convergence process. In order 
to have the system meet the predictions of first order TDPT, v mt should be large enough 
such that a delta function 5(9) is a good approximation of sin 2 h9t/t(ti9) 2 as it appears in 
the derivation of (1301) . Obviously, such criterion depends on the density of the Of about 
9 = 0. We therefore define a lower threshold 

r low := max — \j- (36) 
|0f fc |^Q \0f K \ 

where only 9f k should be considered in (1361) whose corresponding resonant probability tran- 
sitions (i.e. for the case Of = 0) have interaction times of the order of magnitude of r mn 
such that they are relevant for the experiment. In summary, we expect predicted resonances 
to be fully resolved if 

r low ^ r int < r run (37) 

holds. The discussion of concrete driving laws in chapter II I II shows that this criterion is in 
excellent agreement with our numerical simulations. 
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3. Rotating wave approximation 

Interestingly, (|37|) justifies a rotating wave approximation in (124"]) j^ . This allows us to 
calculate approximate population dynamics of the system that, in contrast to (|29|) . conserve 
the total population probability. 

For simplicity, we will assume that there is only one 9f k close to zero. A rotating wave 
approximation simply sets all other terms in ( |24|) that do not contain this 8f- k to zero as 
they are comparatively fast oscillating, such that we are left with with an effective two-level 
system: 



ib n (r) 
ib k (r) 



<--' ■ (uF? k + Df)b k (r) 

coF kn + D kn )b n (r) . 



e -2TTiO? k T 



(38) 



The behavior of such a system is very well understood. A discussion in terms of Bloch 
equations is, for instance, given in 28]. Eqs. (I3"8"j) especially explain why the population 
dynamics in [231 are reminiscent of Rabi oscillations. The effective Rabi-frequency Q c q can 
be calculated (cf. [28]) to be 



tt eS = J (2ir9f) + A\uFp k + Df\ 2 



T: 



(39) 



int. 



which yields a beating period T B of the population dynamics 

rr, 27r 7rr int 
L r := — 



n 



(40) 



In summary, if we assume the system to have initially been in state k, the population 
dynamics of state n is given by: 



Pn{r) 



sin 



2 / 7TT 



1 + {ir9? k Tint) 2 ' 



(41) 



III. NUMERICAL RESULTS AND DISCUSSION 



In this chapter we will present full numerical simulations of driven elliptical billiards and 
analyze the results with the developed perturbation theory of chapter III BL Details to the 
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numerical calculation of the predicted quantities can be found in appendix |A] All numerical 
calculations have been run for r run = 100 periods of driving and h and /i have, w.l.o.g., been 
set to 1. We will always drive the semi-axis a(t) harmonically, i.e. 



a(t) = a + A sin(wt), (42) 

and adjust bit) such that the billiard is driven in different ways as will be specified later. 
To be able to compare the different driving laws, we have chosen to keep the following 
parameters fixed: 



a :-- 



a(t = 0) = 1, A = 0.1 and b := b(t = 0) = V(151. (43) 



These parameters are the same as in 231 ] . 



The energy E(r) will be a key observable for the analysis of the billiard dynamics. It 

2/i 



is calculated from the expectation value of the Hamiltonian H = — f-A as it appears in 



Eq. ©: 

E(t) = ( Mr 



nr)). (44) 

If we apply again the coordinate transformation fl3]) and the unitary transformation U (HJ), 
the energy reads 

E{t) = (A(r) | U\t)H m (t)U(t) \ A(r)> , (45) 

where Hm{t) is given by (JSJ). We can therefore calculate E(r) by determining the population 
p n (r) of the eigenstates of W(t)Hm(t)U(t) in |A(r)), weighting these populations with the 
respective eigenvalues E n (r) of W {t)H m {t)U{t) and sum the results up: 

E(r) = J2E n (r) Pn (r). (46) 

n 

Note that, due to U being a unitary transformation, the energy eigenvalues E n [r) of 
W(t)Hm(t)U(t) are actually identical to the instantaneous eigenvalues E n [r) of H m {t). 

From now on, we refer to the instantaneous eigenstates of U'(t)Hm(t)U(t) as energy 
eigenstates. We understand in particular that the instantaneous eigenstates |n;r(r)) of 
Hm{j) are in general not identical to the energy eigenstates, but unitarily transformed 
energy eigenstates, given by U' \n;r(r)). Note that U is also invariant upon sign-change 
of r) and £, such that an instantaneous eigenstate \n;r(r)) and its corresponding energy 
eigenstate W |n; r(r)) have the same rj- and ^- parity. 
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We will initialize the system in the fourth energy eigenstate (at r = 0) and calculate the 
populations p n (r) upon driving. Note that in [23| instantaneous eigenstates \n;r(r)) were 
used as initial states and for population analyses. 

Interestingly, we find that the overlap | (n; r(r) | U(t) \ n\ r(r)) | 2 of all relevant instan- 
taneous eigenstates |n;r(r)) with their respective energy eigenstates is greater than 94.5% 
for the later analyzed parameter regimes in chapter IIII1 We thus expect that the energy 
eigenstates are similar to the instantaneous eigenstates and also show similar population 
dynamics. Consequently, we will from now on disregard the differences between \n;r(r)) 
and the energy eigenstates W |n; r{r)) when analyzing the billiard dynamics perturbatively 
and will subsequently refer to \n; r(r)) simply as energy eigenstates. This approximation en- 
ables us to predict the seemingly complicated population dynamics of the energy eigenstates 
by our perturbation theory for instantaneous eigenstates. Although the population dynam- 
ics of the energy eigenstates are not expected to be qualitatively different, some predictions 
may be compromised quantitatively. For instance, due to the different actions of U*(t) on 
different |n; r(r)), shifts of the resonance frequencies ( 13Tj) are to be expected. However, as U 
OH) approaches unity for vanishing u, these shifts will rather be observed for larger resonance 
frequencies. We will also find that the resonance shifts become more negligible for higher 
order photon processes. 

In the approximation of instantaneous eigenstates |n;r(r)) being energy eigenstates, a 
transition to a higher excited state increases the energy E(t) (146]) while a transition to a 
lower excited state decreases it. We can therefore determine if a population transition occurs 
at a certain driving frequency u upon simulation time r mn by recording the maximal and 
minimal energy of the billiard in dependence of cu. 



A. Axes-ratio-preserving driving law 

In the following, we provide numerical solutions for various driving laws and analyze them 
with the developed perturbation theory of chapter III Bt A simple but illustrative driving law 
is the so-called axes-ratio-preserving driving law which merely rescales the billiard by varying 
its volume V(r) while keeping the ratio of the semi-axes r(r) constant for all times upon 
driving. As a(t) is given by Eq. ( l4"2j) . we find b(t) = r a(t) and choose r = V0.51 to satisfy 
Eq. (03]). 
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FIG. 1. Energy eigenvalues E n (r) of eigenstates \n;r(r)) with even r]- and £- parity eigenvalues 
for the axes-ratio-preserving driving law. The instantaneous shape of the ellipse at five different 
values of r is drawn below the energy eigenvalue curves. Parameters: ao = 1, &o = \/0.51, A = 0.1. 

The axes-ratio-preserving driving law has the nice property that the instantaneous eigen- 
states |n; r(r)) become time-independent due to fixed r(r) = r , while the eigenvalues E n (r) 
stay time-dependent as can be seen from Eq. f|T3|) . Their variation is solely given by the 
global prefactor 1/V(r) which particularly prevents crossings of energy eigenvalues. Fig. [1] 
shows a sample of eigenvalue curves for one period of driving. The fourth energy eigenstate 
has even rj- and £- parity. Thus, only energy eigenstates in the corresponding sub-Hilbert 
space couple to the chosen initial state. 

In Fig. [2] the dependence of the maximal and minimal energy that has been reached upon 
driving as a solution of the TDSE in Eq. (CEJ) is plotted depending of the driving frequency 
u>. We clearly see sharp peaks and dips at certain driving frequencies. The vertical lines 
represent our predictions of resonance frequencies according to Eq. fl3Tl) . Note that the 
observed resonances deviate slightly from the predicted ones, especially for larger driving 
frequencies. This is due to the unitary transformation U (j4]) that has been neglected in 
our considerations, i.e. we apply perturbation theory only to the instantaneous eigenstate 
\n;r{r)) that is most populated in the energy eigenstate U\t) \n;r{r)). Beside this antic- 
ipated small deviation, we find a very good agreement of the numerical calculations with 
our predictions. All resonances with a comparatively small interaction time r; nt have been 
resolved, while resonances with very large interaction times could not be observed. Natu- 
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FIG. 2. Dependence of the maximal and minimal energy on the driving frequency oj for the axes- 
ratio-preserving driving law. The vertical lines show all predicted resonance frequencies (|3ip with 
an interaction time T mt of less than 2000. The darker the lines are, the longer is the corresponding 
interaction time. Numerical values can be taken from Tab. [J in appendix [Bl Parameters like in 
Fig.ffl 

rally, for interaction times that are longer than (half) the runtime of an experiment r mn , 
a full population transition from the initial state to some other energy eigenstate can not 
happen according to Eq. (jUJ). This is the reason why some peaks in Fig. [2] that correspond 
to transitions to the same energy eigenstate possess different heights. It is interesting to 
note that, although the runtime r run of our numerical simulations was only 100 periods of 
driving, resonances that correspond to an interaction time of up to 2000 periods of driving 
could still be partly resolved in Fig. [2j Tab. [J in appendix [B] provides numerical values 
for all predicted resonance frequencies between and 16 that have an interaction time of 
less than 2000. It also shows that the lower threshold T\ ow is always much smaller than the 
corresponding interaction time r int , such that the first part of Eq. ( 1371) is fulfilled and TDPT 
of first order is applicable. 

One might wonder about a structure of several small peaks, especially for frequencies 
u> > 10.5. We assume that these smaller, not predicted peaks correspond to transitions of 
second order where population is first transferred to one excited state and then from this 
state again transferred to yet another energy eigenstate. This is, for instance, supported 
by a population analysis in Fig. [3] for the small peak at u — 10.81. We see that the 
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FIG. 3. Population dynamics p n i T ) for the axes-ratio-preserving driving law at uj = 10.81. Param- 
eters like in Fig. [TJ 



mean (or envelope behaviour) population of the energetical ground state decreases while the 
population amplitude of the tenth energy eigenstate increases. This may be interpreted as an 
interaction of these two states that consecutively leads to a transfer of population that was 
initially transferred to the ground state and is then pushed to the tenth energy eigenstate. 
Such a process is not included in the time-dependent perturbation theory of first-order in 
chapter III Bl and the dynamics visualized in Fig. |3] are a precursor to the breakdown of this 
simple theory for higher driving frequencies where indirect transitions become more and 
more important. 

It is also interesting that the resonance at uj ~ 15.17 that corresponds to a 4-photon 
transition from the initial state to the 22nd energy eigenstate is so well resolved although 
the interaction time of this resonance is much larger than the interaction time of several res- 
onances that are much worse resolved. The reason for this is that the 22nd energy eigenstate 
has a much higher energy eigenvalue than, for instance, the seventh energy eigenstate. Thus, 
the energy is significantly increased for already a small amount of transferred population 
probability from the initial state to the 22nd energy eigenstate. 
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FIG. 4. Population dynamics p n i T ) for the axes-ratio-preserving driving law at uj = 15.17. Param- 
eters like in Fig. [TJ 

B. Other driving laws 

We have seen in the last chapter that the predictions of TDPT work very well for the axes- 
ratio-preserving driving law. To demonstrate the general applicability of the perturbation 
theory of chapter III we will analyze two further driving laws. The so-called breathing 
driving law b(t) = a(t) — a + bo, where a(t) is again given by Eq. fj42|) . was already discussed 



m 



23| . Fig. [5] shows the eigenvalues of energy eigenstates with even rj- and £- parity for 



one period of driving. We see that crossings of energy eigenvalues are, in contrast to Fig. (TJ 
now possible as r(t) is no longer kept constant. The eigenvalues of the lowest excited states 
are, though, very similar to the ones in Fig. [1] and, consequently, the resonances in Fig. [6] 
resemble the ones in Fig. [2j In Fig. [6] more resonances can be resolved due to a sufficiently 
small interaction time rj nt . This observation can be understood as follows: While for the 
axes-ratio-preserving driving law the transition matrix _D™ m is identical to zero due 
to r(t) = 0, this is not the case for the breathing driving law. The additionally resolved 
resonances for the breathing law correspond, thus, to Landau-Zener transitions with r(t) 



being the Landau-Zener velocity 



321 ] . This role of _D™ m triggering Landau-Zener transitions 



will be even more pronounced for the next driving law that is presented. 

The so-called volume-preserving driving law is just the opposite of the axes-ratio- 
preserving driving law. It keeps the volume V(t) of the billiard fixed, while varying the ratio 
of the semi-axes r(t). Thus, b(t) depends on a{t) (j4"2]) as b{t) = a b /a(t). Fig. [7] shows the 
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FIG. 5. Energy eigenvalues E n {r) of eigenstates |n; r(r)) with even rj- and £- parity for the breathing 
driving law. Parameters like in Fig. [TJ 
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FIG. 6. Analogue of Fig. [2] for the breathing driving law. Numerical values can be taken from Tab. 
HI] in appendix [Bj Parameters like in Fig. [TJ 

corresponding energy eigenvalues. We see that the fourth and seventh energy eigenvalue get 
close upon driving such that we expect that the transition matrix Df m f )23|) couples these 
states strongly. We point out that we can arbitrarily control how close these eigenvalues get 
upon driving by choosing r(t) appropriately. As Df m ( |23l) depends on u only through the 
phase factor e~^ Avmn ^ T \ we expect it to be slowly varying with u, thus setting an upper 
bound on the interaction time (|35j) of resonant population transitions between the fourth 
and seventh energy eigenstate even for small u and corresponding large photon process 
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FIG. 7. Energy eigenvalues E n (r) of eigenstates \n;r(r)) with even r]- and £- parity eigenvalues 
for the volume-preserving driving law. Parameters like in Fig. [TJ 



orders I. 

This expectation is fully confirmed by Fig. [HJ All resolved resonances with u> < 5.5 
correspond exclusively to transitions between the fourth and seventh energy eigenstate. As 
U (S]) gets close to identity for small driving frequencies u, the induced resonance shift 
is negligible. One might wonder, why most of the resonances at small driving frequencies 
have numerically not been fully resolved although their interaction time is short enough. In 
review of Eq. (|32|) . we understand that a detuning in the driving frequency is also multiplied 
by the photon process order I. Thus, one has to adjust the driving frequency very carefully 
to resolve a multiple photon resonance. We assume that the frequency grid in Fig. |8j is not 
fine enough to resolve all predicted resonances although it should in principle be possible. 

Finally, population analyses close to the resonance frequencies confirm the Rabi-like be- 
havior of the population dynamics as predicted by Eq. 04 II) . This can be especially well 
illustrated for resonances with small interaction times. As an example, Fig. [9] shows almost 
perfect Rabi-like population dynamics of the fourth and seventh energy eigenstate in the 
nearly resonant case of a three photon process at u = 3.32. Comparison of the observed 
beating periods with the corresponding interaction times r int in Tab. IIHI in appendix [B] 
gives, even quantitatively, a very good agreement as predicted by Eq. (JlDjl . We point out 
that as we can tune the strength of the transition matrix _D" m by choosing how close the 
energy eigenvalues get upon driving, we can also tune the interaction time r int in the regime 
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FIG. 8. Analogue of Fig. [2] for the volume-preserving driving law. Numerical values can be taken 
from Tab. IIIII in appendix [Bj Parameters like in Fig. [TJ 
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FIG. 9. Population dynamics p n (T) for the volume-preserving driving law at uj = 3.32. Parameters 
like in Fig. [TJ 

of weak driving where it is mainly determined by D" m . Hence, we can, in principle, also 
control the beating period of the present effective two-level Rabi-system. 



IV. BRIEF SUMMARY 



A time-dependent perturbative approach for elliptical quantum billiards with oscillating 
boundaries has been developed. As our major results we have obtained a Fermi Golden 
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Rule, predicting the driving frequencies yielding resonant population transfer between in- 
stantaneous eigenstates as observed in [23j and a criterion allowing to decide which of these 
resonances are observable in a corresponding experiment of certain duration. Extensive 
numerical simulations have been performed for three different driving laws which are in ex- 
cellent agreement with our predictions. Particularly for the volume preserving driving law, 
due to the change of the billiard geometry upon driving, Landau-Zener transitions have been 
analyzed to take place. Depending only weakly on the driving frequency, these transitions 
allow for resonant population transfer also for very weak driving. We have shown, that 
the billiard dynamics can be reduced in this regime to an effective two-level system with in 
principle arbitrarily tunable oscillation period. 

Further interesting phenomena beyond the scope of our perturbative description can be 
expected in the numerically challanging regime of strong driving. 
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Appendix A: Matrix elements 



Introducing the matrices 



/'= E 

n,m,n', 

P= E 

n,m,n' , 

P= E 

n,m,n', 

/*= E 



n,m,n' ,m' 



f= E 

n,m,n',r 

/•= E 



^n,m) ^m,m' fnmn' (^n',r. 



^n,m) $m,m' fnmn' {^n',r. 



&n,m) 3(m—2),m' fnmn' i^n' 
$n,m) 5{m-2),m' fnmn' (®n' 
^n,m 

>m' fnmn' ^^ti' 
&n,m) $ (m+2) ,m' f nmn' i^n' 



n,m,n' ,m' 



with matrix elements 



f 1 

J nmn 
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J n 
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n m,n r 

^ On,n' 



f2 

' nmn' 



3 

nmn' 



^Jm+1 y^m,n) Jm+1 (^m,n') Jo 



"^71-2,?!.' 
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(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



(A7) 
(A8) 

(A9) 
(A10) 

(All) 
(A12) 



4</ m +l(^m,n)^m+3(^m+2,n') Jo 

where J m is again the cylindrical Bessel function of order m and fc min is its n-th root. We 
have a convenient form of representing H u (ED , M(r) ( fT2l) and Hf(t) ( jTUi) in the eigenbasis 
{|$n,m)}n,m (13) of the static circular billiard: 



HM = 9i(T)f 1 + g 3 (T)(f s + f 



M(r) 



(A13) 



(A14) 
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H f {t) = g 2 {r)f 2 + g A (r) (/ 4 + f) . (A15) 

Diagonalizing M(r) yields the instantaneous eigenstates \n;r) and their eigenvalues q n {r). 
One could in principle calculate the energy eigenvalues E n {r) = M ^ T - ) g n (r(r)) from the 
g n (r), but it turns out that diagonalizing Hm directly increases the numerical precision of 
the energy eigenvalues. Note that the time-dependent factors gi(r) as well as the matrix 



elements fl mn i are the same as in 



231 ] . However, the matrix elements f nmn > have been 



reduced to a much simpler form, using orthonormality relations of the Bessel functions. 



h 2 



01 T 



H \cl(t) 2 b(r) 2 
92(r) = y. (a(T)a(r) + 6(r)6(r) 
1 1 



h 2 
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9A T ) 



fi \a(r) 2 b(r) 2 
jj, (a(r)a(T) — b(r)b(r 



(A16) 
(A17) 
(A18) 
(A19) 



Further note that the sign of (r) is inverted in comparison with [23]. We can now calculate 
the transition matrix elements _D" m (123!) and F" m (|22j) : 
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jpnm 
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u i ' u 2,l 
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; 4J 
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nm 
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After calculating these quantities and diagonalizing M(r) and Hm, it is straightforward to 
reproduce the theoretical predictions contained in this work. 

Appendix B: Theoretical predictions 
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TABLE I. Numerical values for all predicted resonance frequencies between the values and 16 
with an interaction time T m t less than 2000. Information is provided on the corresponding lower 
threshold T\ ow , the quantum number of the coupling instantaneous eigenstate n and the photon 
process order I of the resonance. 
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TABLE II. Analogue to Tab. Ubut for the breathing driving law. 
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